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Lattice discretization of the supersymmetric Yang-Mills quantum mechanics is reviewed and 
results of the Monte Carlo simulations of the simplified model are presented. The D=4, N=2, 
quenched system, studied at finite temperature, reveals existence of the two distinct regions 
which may correspond to a black hole and the elementary 0-branes phases of the M-theory 
conjectured in the literature. New results for higher gauge groups A'' < 9 lead to the similar 
picture, however the nature of the transition between the two phases and its precise scaling 
with A'^ is yet unresolved. 

1 Matrix quantum mechanics and its lattice formulation 

Recent developments in nonperturbative string theories have lead to an exciting hypothesis of 
the existence of a single theory (M-theory) which encompasses five known superstrinff theories 
and eleven dimensional supergravity in a unified scheme (for a recent review see e.g. u). Even 
though its details are not known, the M-theory has a tantalizing potential to unify all interactions 
and particles. In particular it may offer a topological explanation of such fundamental features 
as three families and fractional charges. It may lead to a standard model gauge group with 
AA=1 supersymmetry. It provides understanding of the Bekenstein-Hawking entropy puzzle and 
much more. Moreover, Banks, Fishier, Susskind and ShenkerB suggested that the spectrum of 
M-theory is equivalent to that of a supersymmetric Yang- Mills quantum mechanics (SYMQM) 
which results from the dimensional reduction of the 10 dimensional supersymmetric Yang-Mills 
theory. This allows to use a host of nonperturbative methods to quantitatively study both 
systems. Accordingly we have constructed the Wilson discretization of the above quantum 
mechanics and proposed to investigate it with the standard lattice techniques u. The ultimate 
goal is to study D=10 SYMQM for the large size of the SU(N) matrices. However, even for 
this relatively simple, one dimensional quantum mechanical system this is a formidable task, 
the main difficulty being the complex fermionic determinant (and pfaffian) in D=10, and time 



consuming procedures for simulating nonabelian systems of large matrices. On the other hand 
the system can be simulated at present in other regions of the parameter space {D, N, Nf) and 
such a study may provide some information about its general structure. Therefore we have 
decided to set up a systematic lattice study of SYMQM beginning with the simplest case of 
D = 4,N = 2, Nf = and gradually extending it as far as possible towards the BFSS limit. 
In fact in this talk I will also report on some results for higher N. Supersymmetric Yang-Mills 
Quantum Mechanics has been intensively studied a. Although its exact solution is still not 
available, many results are known and a) can be tested and b) provide us a guidance from a 
"simple corner" of parameter space to the ultimate "BFSS corner". 
The action of the SYMQM reads 

S = Jdt QT¥F^,(f)2 + ^''{t)V^''it)^ . (1) 

where /i, = 1 . . . D, all fields are independent of the space coordinates x and supersymmetric 
fermionic partners belong to the addjoint representation of SU(N). The discretized system is 
put on a D dimensional hypercubic lattice A'^i x . . . x Nd reduced in all space directions to 
Ni = 1, i = 1 . . . D — 1. Gauge and fermionic variables are assigned to links and sites of the new 
elongated lattice in the standard manner. The gauge part of the action reads 

m . 

Sg = -PT.T. n^^"^ U^uim)), (2) 

m=l fJ.>i> 

with P = 2N/a^g'^, and U^u{m) = Ul{m)Uj^{m+v)Uv{m+^)U^{rn), U^{m) = exp {iagA^{am)), 
where a denotes the lattice constant and g is the gauge coupling in one dimension. The integer 
time coordinate along the lattice is m. Periodic boundary conditions U^{m + v) = U^{m), 
u = 1 . . . D — 1, guarantee that Wilson plaquettes U^u tend, in the classical continuum limit, to 
the appropriate components F^^, with space derivatives absent. In this formulation the projection 
on gauge invariant states is naturally implemented. Eq. (§) is the basis of the MC simulations of 
the still more simplified, D=4, N=2 and Nf = (quenched), model Some preliminary results 
for higher gauge groups < 9 are now also available and will be shortly discussed. 

2 Results 

As a first problem we have chosen the question of the phase structure of the model. A very 
interesting feature of the M-theorv is the solution of the Bekenstein-Hawking entropy puzzle 
in terms of the elementary branes a. In particular, the theory predicts existence of the phase 
transition between a low temperature "black hole" phase and a high temperature phase described 
by the elementary branes Q ^. Threfore a natural question arises if the simplified model (i.e. 
SYMQM) posesses any nontrivial phase structurei. At the same time it is well established that 
QCD (or pure Yang-Mills theory) has two different phases (e.g. confinement and deconfinement). 
Since the action (|2|) is basically QCD-like, one may expect that the dimensionally reduced model 
may indeed exhibit similar phenomenon. To check this we have measured the distribution of 
the trace of the Polyakov line 

^ = ^Tr[n (3) 

\m=l / 

which is a very sensitive determinant of the two phases in QCD. Similarly to lattice pure gauge 
system, symmetric concentration of the trace around zero indicates the low temperature phase 
with < P >= (here a "black hole" phase) while clustering around ±1 (or, for arbitrary N, 

° The full phase structure of the M-theory is expected to be much more complex op. cit. 



around the elements of Z^) is characteristic of the high temperature phase (here the elementary 
0-brane phase). 

There is one important difference between the dimensionally reduced model and the pure 
gauge system in the large spatial volume. Of course the one dimensional system with local 
interactions cannot have a phase transition for finite N. However at infinite N the sharp singu- 
larity may occur. This suggests that the number of colours plays a role of a volume (indeed the 
number of degrees of freedom is proportional to A^^). Consequently, and similarly to a small 
volume systems, we expect to see some signatures of a phase change for finite and even small 
N . However the generic singularities would develop only at infinite N. In fact studying the N 
dependence of the above signatures allows often to determine the properties of the infinite N 
phase transition via the final size scaling analysis. 

In ReftI above program was carried out for D=4, N=2 and in the quenched approximation 
Nf = 0. It was found that the distribution of (^ is indeed changing from a convex to a concave 
one at some finite value of /3 = Pd^t)- Moreover the dependence of Pc on Nt turned out to be 
consistent with the canonical scaling. The best fit to Monte Carlo results gave Pd^t) = (0.17 it 

f 3 02^z0 33^ "i 

0.05)NI: ' in good agreement with the expected dependence Pc ~ ^t- This indicated 

that in the continuum limit the transition occurs at finite temperature Tc = A{g'^N)^^^^\ with 
A = .28 lb 0.03. We have also measured the average size of the system = X]a(^?)^ which is 
shown in Fig.l for different values of the temperature. Our results agree qualitatively with the 
mean field calculation at large Nu. Moreover, the pseudocritical temperature determined from 
Fig.l is consistent with the value of A quoted above, which to our knowledge was determined 
for the first time. Obviously there is a long way between N = 2 and the BFSS limit N — > oo. 
Even though there is an intriguing evidence from lattice simulations that large N limit may be 
reached rather early in N (at least for some observables) iH, it is very important to repeat this 
calculation for higher N. In particular the coefficient A may depend on A^. 



Figure 1. Extrapolated to the con- 
tinuum limit size of the system as a 
function of the temperature (in the 
units of the one dimensional gauge 
coupling g). 
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Currently we are extending this analysis to higher gauge groups (and still for Nf = 0). 
Detailed account of our results will appear elsewhere. Main points are the following. For N > 2 
the trace P in Eq.(^ is complex. The condition det(U)=l constraints the trace to lie inside 
the "N-star" in the complex plane (see Fig. 2). In the low temperature phase the distribution is 
symmetric and peaked around P = ensuring < P >= 0, while in the high temperature phase 
one expects P's to be concentrated around the elements of Z]\f. A sample of our preliminary 
results for N=3,5,8 and at high and low temperatures (/? ~ T^) is shown in Fig. 2. Indeed 
one clearly sees existence of both (high and low temperature) regimes. Moreover, in the high 
temperature region, the system has a tendency to be stuck in one of the Zf^f minima. This is a 
typical indication of the spontaneous symmetry breaking which might occur at ^ oo leading 
eventually to a non-zero value of the order parameter < P >. We expect to accumulate soon 
enough data to be able to decide if the two above regions are separated by the genuine phase 
transition or a smooth cross-over. 

To conclude, the lattice simulations provide a novel approach to the Yang-Mills quantum 
mechanics and possibly to the M-theory. Many problems (e.g. identification of the graviton 
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Figure 1: Distribution of the Polyakov line in the low (upper row) and high (lower row) temperature phases of 
the quenched Yang-Mills quantum mechanics with = 3, 5, 8. 

multiplet) can be formulated and quantitatively attacked in this framework. On the other hand 
much more work is required to continue development of this fascinating subject. In partic- 
ular, construction of the new algorithms capable to deal effectively with dynamical fermions 
propagating in this one dimensional system is an open question. 
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